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Differential cross sections of the Coulomb deexcitation in the collisions of excited muonic hydrogen 
with the hydrogen atom have been studied for the first time. In the framework of the fully quantum- 
mechanical close-coupling approach both the differential cross sections for the nl —> n'l' transitions 
and Z-averaged differential cross sections have been calculated for exotic atom in the initial states 
with the principle quantum number n = 2 — 6 at relative motion energies i?cm = 0.01 — 15 eV and 
at scattering angles 6cm = 0 — 180°. The vacuum polarization shifts of the ns-states are taken into 
account. The calculated in the same approach differential cross sections of the elastic and Stark 
scattering are also presented. The main features of the calculated differential cross sections are 
discussed and a strong anisotropy of cross sections for the Coulomb deexcitation is predicted. 


Introduction. Exotic hydrogen-like atoms are formed 
in excited states, when heavy negative particles {p~ , tt~ , 
etc.) are slowed down and captured in hydrogen media. 
The following atomic cascade of collisional and radiative 
transitions proceeds via many intermediate states up to 
nuclear absorption or transition to the ground state oc¬ 
curs. Since the experimental data are mainly available for 
the last stage of this atomic cascade, the reliable knowl¬ 
edge of the total and differential cross sections (DCS) of 
the collisional processes during the cascade is needed for 
the realistic analysis of these data. 

In particular, the collisional processes 

{p~p)nl + His {p~p)n'l' + His ( 1 ) 

of the elastic scattering (n' = n, V = 1), Stark transitions 
(n' = n, and Coulomb deexcitation {n' <n) essen¬ 

tially change the energy- and nCdistributions of exotic 
atoms. It is especially important from the view of the 
precise experiments at the Paul Scherrer Institute with 
muonic Q and pionic 0 hydrogen atoms. These ex¬ 
periments are aimed at the extraction of the root mean 
squared proton charge radius with the relative accuracy 
of 10“^ from the Lamb shift experiment Q and the deter¬ 
mination of the ttN scattering lengths with the accuracy 
better than 1% by extracting the shift and width of Is 
state due to the strong interaction in pionic hydrogen Q . 
The proper analysis of these experiments requires a reli¬ 
able theoretical treatment of related cascade processes. 

The Coulomb deexcitation (CD) plays an important 
role in the kinetic energy distribution of the exotic atoms. 
In particular, the energy distribution of {pTp)- and 
[t:~p)- atoms during the radiative transitions np —> Is 
have significant high-energy components resulting from 
the preceding Coulomb deexcitation transitions. Before 
the recent paper Q the only An = 1 transitions are as¬ 
sumed to be important in the Coulomb transitions at 
low n. Moreover, it is also suggested in the cascade cal¬ 
culations 0 that the CD process results in the isotropic 
angular distribution. 

The main goal of this paper is to introduce the first 
theoretical study of the Coulomb deexcitation differen¬ 
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tial cross sections in the muonic atom - hydrogen atom 
collisions. In particular, we are interested in the main 
features of the n and E dependences of these cross sec¬ 
tions. Our new results, concerning the DCS of elastic 
and Stark scattering are also presented. 

The first theoretical study of DCS for the elastic and 
Stark scattering in the collisions of the excited exotic 
atom from atomic hydrogen has been performed within 
the quantum-mechanical adiabatic approach [^. Later, 
these cross sections were also calculated in the framework 
of the close-coupling model and semiclassical approxima¬ 
tion where the Coulomb interaction of the exotic atom 
with the hydrogen atom field is modeled by the screen¬ 
ing dipole approximation. This approximation becomes 
invalid at low collisional energies below some value E*, 
which depends on the principal quantum number n (see 
0 )- Thus, the latter approach as well as various mod¬ 
ifications of the semiclassical model can result in 

uncontrolled errors in the low-energy region where only 
a few partial waves are important. 

In the present paper the DCS of the processes (1) are 
studied in the framework of the more accurate close cou¬ 
pling ( CC) approach. The approach has been developed 
earlier by the authors to describe the elastic scatter¬ 
ing and Stark transitions in the exotic atom - hydrogen 
molecule collisions and employed recently 0 to give a 
unified treatment of the processes (1). In the framework 
of the CC approach the first fully quantum-mechanical 
calculations of the total cross sections of the Coulomb de¬ 
excitation both for muonic 0 and pionic m hydrogen 
atoms have been performed. 

Approach. In the framework of the close-coupling 
approach the total wave function of the four-body sys¬ 
tem is expanded in terms of the basis states with the 
conserving quantum numbers of the total angular mo¬ 
mentum JM and parity tt = (—1)*+^. The basis states 
are chosen as tensor products of the corresponding wave 
functions of the free exotic and hydrogen atoms and the 
angular wave function of their relative motion. This ex¬ 
pansion results in the set of the coupled differential equa¬ 
tions. In contrast with paper in the present approach 
the interaction potential matrix for the exact four-body 
Coulomb interaction of the colliding atoms is calculated 
analytically. Moreover, as it is shown in paper Q the 
approximation similar to ’’dipole approximation” results 
in an improper description of the CD process, especially 
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for the low-lying states of the exotic atom. 

At fixed Ecm (energy of collision in the center of mass 
system) and given J and tt the set of the coupled equa¬ 
tions are solved numerically by the Numerov method 
with the standing-wave boundary conditions involving 
the real and symmetrical AT-matrix. The corresponding 
T-matrix can be obtained from the K matrix using the 
matrix equation T = 2iK (/ — i K) . In the present 
study, as distinct from paper Q, we take into account 
the energy shifts of the ns states due to the vacuum po¬ 
larization which are very important especially at the low 
kinetic energies comparable with the energy shift and for 
the lower states of the exotic atom. The 2s — 2p energy 
splitting Ae 2 s- 2 p for the muonic hydrogen atom is equal 
to 0.206 eV and this splitting decreases approximately as 
with n increasing. 

The formalism has been described in more details 
in 0. Here we give only the formulas for the DCS. The 
differential cross sections for the transition from the ini¬ 
tial state {nl) to the final state {n'V) are defined as 


d(Jnl- 


dVt 


1 fc/ 

21 + 1 ki 


^ ( I/nZm—(^2; H) I , (2) 


where the scattering amplitude for the transition nlm 
n'Vm' is given by 
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states with I ^ 0 in the entrance channel, which are as¬ 
sumed to be degenerated. 

Results. The numerical calculations of the DCS 
for the collisional processes (1) have been done for 
(/ip)„ atoms with the initial principal quantum num¬ 
ber values n = 2 — 6 at the relative motion energies 

Ecm = 0.01 — 15 eV and at all the scattering angles 
0cm from zero up to 180°. All the exotic-atom states cor¬ 
responding to the open channels have been included in 
the close-coupling calculations. Some of our results are 
present here in Figs. 1-9 both for the cross sections of 
the separate nl —> n'V transitions and for the Z-averaged 
ones. 



FIG. 1: Differential 21 21' and 21 —> Is cross sections for 

(/rp) 2 i + H collisions vs. cms scattering angle 9cm at Ecm = 
0.01 eV (referring to 2p state). 


Here, ki and kf are the cms relative momenta in the ini¬ 
tial and hnal channels, correspondingly; VI = 9^(p, where 
9 is the cms scattering angle, and is the transition 

matrix in the total angular momentum representation. 
The indices of the entrance channel and the target elec¬ 
tron state are omitted for brevity. 

In order to illustrate the most general features of DCS, 
it is also useful to introduce the cross sections averaged 
over the initial orbital angular momentum I of the exotic 
atom. So, the following Z-averaged angular distributions 
are also discussed in the present study: 
for elastic scattering 




^(2Z -I-1) 


dVt 


(4) 


for Stark transitions 



dfl 




(5) 


and, finally, for CD process 



dfl 


/ 1 ' 


n < n. (6) 


Hereafter, the atomic units will be used throughout the 
paper and the collision energy will be referred to the 


In Figs. 1 and 2 the DCS for 2p —> 2p, 2s (elastic and 
Stark scattering) and 2p Is (Coulomb deexcitation) 
transitions versus cms scattering angle 0cm at energies 
Ecm = 0.01 eV and 1 eV (referring to the 2p threshold) 
are shown. It is worthwhile noting that the relative mo¬ 
tion energy in the entrance channel for the 2s ^ 2s, 2p 
and 2s —> Is transitions increases due to the Lamb shift 
A£ 2 s- 2 p = 0.206 eV in comparison with the scattering 
processes of the muonic hydrogen atom in the 2p state. 
As it is seen from Fig. 1, the angular distributions of 
the elastic 2p —> 2p scattering and Coulomb 2p ^ Is 
transition are similar and almost isotropic, their shapes 
are mainly defined by the contributions of the S-wave 
relative motion with a small mixture of the P-wave. 

It is also seen that Coulomb deexcitation process for 
2p —> Is transition is more than four order of the magni¬ 
tude suppressed relatively to the 2p —s- 2p elastic scatter¬ 
ing one and more or about two order of the magnitude 
as compared with both the Stark 2p 2s and Coulomb 
2s —> Is transitions, respectively. 

The Stark transitions 2p —> 2s and 2s —> 2p are 
strongly inelastic at such a low energy (the initial and 
final energies differ more than twenty times) and that 
results in their quite unusual angular distributions (see 
Fig. 1) with the maximum at zero scattering angle and 
minimum at 0cm ~ 30° and smooth increasing at the 
backward hemisphere. 

In the DCS of the elastic 2s —> 2s scattering (Acm = 
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FIG. 2: The same as in Fig. 1 but at Ecm = 1 eV. 


0.216 eV) we see the structures (due to higher partial 
waves, involved in the scattering process) as compared 
with the almost isotropic DCS for 2p 2p transi¬ 
tion. The comparison of the DCS of the Stark transi¬ 
tions, presented in Figs. 1 and 2, shows a significant 
effect of the Lamb shift on these cross sections, which 
changes substantially depending on the ratio value of the 
i?cm/Ae 2 s- 2 p- The Stark DCS of the 2p^2s and 2s^2p 
transitions at low kinetic energy comparable with the 
Lamb shift value have quite different behavior in compar¬ 
ison with the elastic ones and reveal a strong suppression 
in the forward hemisphere by more than three order of 
the magnitude. When the energy of the collision is much 
larger than Ae 2 s- 2 p we observe the usual picture iH 
of the DCS for Stark transitions with a strong forward 
peak and a set of maxima and minima (see Fig. 2). 

Coulomb deexcitation process occurs at essentially 
smaller distances than the elastic and Stark processes. 
So, the number of partial waves involved in the CD pro¬ 
cess, as a rule, is much smaller, as a strong centrifugal 
barrier prevents the colliding atoms from penetrating in 
the interaction region corresponding to the process. With 
decreasing the value n the number of the partial waves 
contributing to the CD process (at the fixed energy) also 
decreases. 

In particular, the DCS of the CD process in case of 
n = 2 have a quite simple angular dependence due to a 
few partial waves involved in the process( see Fig. 1 and 
2). The angular dependences of the DCS for 2s —> Is 
and 2p —>■ Is Coulomb transitions are mainly deter¬ 
mined by the contributions of the lowest partial waves 
of the relative motion at all energies under considera¬ 
tion. The differential cross sections of these Coulomb 
transitions as it is seen in Fig.2 have a similar angular 
dependence which shape is slowly changed enhancing the 
forward hemisphere scattering with energy increasing. 

It is well-known that in the muonic hydrogen atoms, 
the 2s-state plays a particular role due to 2s Lamb shift 
and has no analog in the other exotic atoms in which the 
strong interaction leads to a large rate of the nuclear ab¬ 
sorption from this state. In particular, a new knowledge 
about the collisional quenching of 2s state at collisional 
energy near or below the 2p threshold is of special inter¬ 


est. 

As it shown in Fig. 1, the DCS for the 2p—>ls transi¬ 
tion is strongly suppressed in comparison with the main 
2s —> Is transition about two order of the magnitude 
and this suppression is also observed at higher kinetic 
energy (see also Fig. 2). Hence the 2s ^ Is transition 
determines the CD 2^1 transition at all kinetic ener¬ 
gies and it is quite probable that the observed collisional 
quenching of the metastable 2s state of the muonic atom 
and the high energy component of muonic hydrogen in 
Is state can be explained by the direct Coulomb deexci¬ 
tation process. 



FIG. 3: Differential elastic {I — I') cross sections for {p,p)u+H 
collisions vs. cms scattering angle 0cm at Acm = 1 eV. 

The typical angular distributions for the individual 
nl —> nV transitions for n = 5 are shown in Fig. 3 
and Fig. 4 for the elastic scattering and Stark transi¬ 
tions, respectively. It is well known that DCS of 

these processes are similar to the diffraction scattering 
(at the collisional energies more or about 1 eV) with a 
strong forward peak which is enhanced with increasing 
energy and a set of maxima and minima. While the elas¬ 
tic cross sections always have a strong peak at 0cm = 0, 
the first maximum position in the Stark DCS depends 
on the = 1^ — l'\ value. In particular, for AZ = 1 this 
maximum is at finite scattering angles as it also remarked 
in [^. According to our calculations, the sharpest vari¬ 
ations in DCS are always observed in the ns —> n's and 
ns n'p transitions (see also Fig. 7 for DCS of the CD 
process). 

The first and next peaks in the forward hemisphere for 
the elastic scattering (see Fig. 3) and Stark transitions 
(see Fig. 4) have a tendency to be less pronounced and the 
angular distribution becomes smoother with increasing I 
and Al is increased. According to our calculations the 
shape of the peaks in the forward hemisphere is sharper 
with n increasing at the fixed collisional energy. 

The dependence of the Z-averaged DCS for the elastic 
scattering on the collisional energy is shown in Fig. 5 for 
n = 4. While the DCS for the individual elastic nl nl 
transitions (the same is valid for the Stark nl —> nV tran¬ 
sitions) reveal the complicated structure, the Z-averaged 
cross sections smooth out many details and allow to study 
the most general features of the process. Figure 5 shows. 
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FIG. 4: Differential Stark cross sections for 5/ ^ 5/' transi¬ 
tions with AZ = 1 — 3 in (^p) + H collisions at Eon = 1 eV. 



FIG. 5: The /-average differential elastic cross sections for 
(/ip)n =4 + H collisions at different energies. 


that at low energies (< 1 eV) the DCS can be approxi¬ 
mately considered as constant for the simple estimation 
at scattering angles 9cm more or about 75°. However, 
at higher energies the appreciable enhancing of the back¬ 
ward scattering is observed. 

In Fig. 6 we show the /-averaged Stark DCS at Ecm = 
1 eV for different values of n. Here one can see that 
with the increase of n the first forward peak also becomes 
sharper and narrower (as well as for the elastic scattering) 
but remains to be at the finite values of scattering angle. 
The height of this peak depends on n not so strong as the 
diffraction maximum in elastic scattering (cf. Figs. 5 and 
6). The diffraction structure of minima and maxima be¬ 
comes less pronounced with increasing n. As a whole our 
results for the elastic and Stark DCS are in a qualitative 
agreement with the previous calculations @ . 

Now we are coming to the discussion of the typical 
angular distributions for the CD process. As far as we 
know, the calculations of these DCS have not been re¬ 
ported until now and in the cascade calculations the an¬ 
gular distributions of the CD process are presumed to be 
isotropic. The calculated DCS for individual nl —> n'V 
transitions with An = 1 and 2 at relative motion energy 



FIG. 6: The /-average differential Stark cross sections for 
{pp)n + H collisions (n = 2 — 6) at Ecm = 1 eV. 


A'cm = 1 eV are shown in Figs. 7 and 8, respectively. In 
Fig. 9 the /-averaged DCS for the 6^5 transition at 
different values of the relative motion energy from 0.01 
up to 15 eV are presented. 

Our study reveals the following main features of the 
CD angular distributions. The angular distributions 
both of the individual and /-averaged cross sections (ex¬ 
cluding very low energies) are far from isotropic: as a 
whole the scattering at 9cm ^ 60° and 9cm > 120° is 
noticeably enhanced. 



FIG. 7: Differential CD cross sections for the individual tran¬ 
sitions with An = 1 for n = 5 at Ecm = 1 eV. 

The DCS for ns —> n's transitions (see Figs. 7, 8) 
have (as in case of the elastic scattering) a more pro¬ 
nounced diffraction structure with sharp maxima and 
minima and a strong peak at zero angle as compared 
with the smoother angular dependence for the other CD 
transitions. This behaviour can be simply explained by 
the conditions L = L' = J (for ns —> n's transitions) 
which strongly reduce the number of terms in the ampli¬ 
tude (3) in contrast with the other transitions. 

The increase of kinetic energy enhances asymmetry in 
the angular dependence of the /-averaged DCS and de¬ 
creases the role of the backward scattering (see Fig. 9). 

Summary. The fully quantum-mechanical CC ap- 
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FIG. 8: The same as in Fig. 7 but for the An = 2 transitions. 



FIG. 9: The /-average CD differential cross sections for tran¬ 
sition 6 —» 5 at different energies. 


proach has been applied for the calculations of the elas¬ 
tic scattering, Stark transition and Coulomb deexcitation 
DCS in a self-consistent manner and the detailed analysis 
of the obtained results has been performed. For the first 
time the DCS of the CD process have been calculated 
for the values of the principal quantum number and ki¬ 
netic energy relevant for kinetics of the atomic cascade. 
The first results for the direct collisional quenching of 
the 2s-state due to CD process were also obtained. The 
present study reveals the new knowledge about the CD 
process and is very important for the reliable analysis of 
the K X-ray yields and high energy component in the 
kinetic energy distribution of muonic hydrogen atoms. 
We hope that our study allows to remove some uncer¬ 
tainties inherent in the previous cascade calculations, 
which resulted from the treatment of collisions, especially 
Coulomb deexcitation, involving different and not always 
self-consistent approximations. 
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